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MODELING AND SIMULATION OF OPTICAL 
CHARACTERISTICS IN A TEXTURED A-SI THIN FILM 
SOLAR CELL USING THE TRANSFER MATRIX METHOD 
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Abstract. A simulation of a thin film hydrogenated amorphous silicon cell was 
developed, in order to demonstrate the use of the transfer matrix method, TMM , with 
diffusive interfaces to model a textured a-Si solar cell. A quantum efficiency comparison 
between experiment and simulation of thin film a-Si:H solar cell was performed. We have 
implemented a self-consistent optoelectronic model for simulation of solar cells with 
rough textured interfaces. 
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1. Introduction 


Worldwide growth of photovoltaics has been fitting an exponential curve for more 
than two decades [1-3]. During this period of time, photovoltaics (PV), also known 
as solar PV, has evolved from a pure niche market of small scale applications 
towards becoming a mainstream electricity source [4]. There that is being made a 
rapid progress with inorganic thin-film photovoltaic (PV) technologies, both in the 
laboratory and in industry. While amorphous silicon based PV modules have been 
around for almost 30 years, recent industrial developments include the first 
polycrystalline silicon thin-film solar cells on glass and the first tandem solar cells 
based on stacks of amorphous and microcrystalline silicon films [5-7]. Optical 
modeling has become a powerful tool for analyzing the optical properties of the 
thin-film solar cells. To model the complex optical behavior in the multi-layer 
optical systems with rough interfaces, such as a-Si solar cells, accurately, there 
should be taken into account both the interference effects and light scattering [8]. 
To optimize the performance of the solar cells by taking advantage of the enhanced 
light absorption, one needs to understand the influence of light scattering on the 
absorption profile in the solar cell. We would expect that the light is completely 
scattered, or in other words, diffused in the solar cell after passing through several 
rough interfaces. However, in wavelength-dependent spectral response and total 
reflection measurements on the solar cells, a moderately pronounced interference 
pattern can be observed. This indicates that besides the scattered light, which 
propagates incoherently in different directions, there is still a significant amount of 
specular light in the structure that propagates coherently. 
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To model the light scattering process in the solar cells accurately, the optical model 
needs to be calibrated well, i.e., the input scattering parameters of the model have to 
reproduce the available experimental data as precise as possible [9]. 


2. Materials and method 
2.1. Materials 


Amorphous silicon (a-Si) has been used as a photovoltaic solar cell material for 
devices which require very little power, such as pocket calculators, because their 
lower performance compared to conventional crystalline silicon (c-Si) solar cells is 
more than offset by their simplified and lower cost of deposition onto a substrate. 
The first solar powered calculators were already available in the late 1970° [10]. 


More recently, improvements in a-Si construction techniques have made them more 
attractive for large-area solar cell use as well. Here their lower inherent efficiency is 
made up, at least partially, by their thinness — higher efficiencies can be reached by 
stacking several thin-film cells on top of each other, each one tuned to work well at 
a specific frequency of light. This approach is not applicable to c-Si cells, which are 
thick as a result of their construction technique and are therefore largely opaque, 
blocking light from reaching other layers in a stack. 


The main advantage of a-Si in large scale production is not efficiency, but cost. 
Amorphous silicon cells use only a fraction of the silicon needed for typical c-Si 
cells, and the cost of the silicon has historically been a significant contributor to cell 
cost. However, the higher costs of manufacture due to the multi-layer construction 
have, to date, made a-Si unattractive except in roles where their thinness or 
flexibility are an advantage. 


Glass Substrate 


Fig. 1. Schematic of a state-of-art p-i-n a-Si:H solar cell on a glass superstrate [12]. 
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Typically, amorphous silicon thin-film cells use ap-i-n structure. Typical panel 
structure includes front side glass, TCO, thin-film silicon, back contact, polyvinyl 
butyral (PVB) and back side glass. Uni-Solar, a division of Energy Conversion 
Devices produced a version of flexible backings, used in roll-on roofing products. 


The structure of a state-of-art p-i-n a-Si:H solar cell on a soda-lime glass superstrate 
is shown in Figure 1. The front TCO (ZnO:Al) must have a low sheet resistance of 
about 10 Q/sq and a high optical transmission for visible wavelengths [11,12]. 


2.2. Method 


Texturing is often used to improve light trapping in thin film solar cells. [13,19]. 
The textured surfaces act to scatter the light in oblique angles lengthening the 
optical path and effectively trapping a greater percentage of light within the thin 
layers [14,15]. 


Numerical simulation including self-consistent inclusion of optical and electrical 
simulation has been an important tool in optimization of solar cell designs. 
However texturing is problematic for numerical optoelectronic simulation of thin 
film solar cells. Purely geometric simulation such as Ray Tracing Method (RTM) is 
not suitable since it ignores interference effects occurring in devices with material 
layers with thicknesses on the order of the wavelength of the incident light. Also it 
is difficult to properly describe randomly textured surfaces geometrically [16, 17]. 
On the other, hand coherent modeling using Transfer Matrix Method (TMM) 
ignores the effects of surface texturing. Finite difference time methods could also be 
used but at an enormous sacrifice in computation time and with problems with 
periodicity and non-normal incidence. 
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Fig. 2. Relationship between the haze functions and the specular and diffused light at a textured 
interface. 
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An approach that accounts for both the coherent propagation of specular light as 
well as the diffuse propagation of incoherent light occurring at the rough interfaces 
in the device has been proposed [8, 18]. In this approach, coherent light is modeled 
using the transfer matrix method while incoherent light is modeled using geometric 
methods. This model is incorporated into the Luminous simulator in the Silvaco 
ATLAS framework [7] for self-consistent optical propagation and drift-diffusion 
simulation of photo-detection problems. 


The relation between coherent and incoherent light is defined by a set of haze 
functions defined at each textured interface. The haze functions Hy and Hp as 
illustrated in Figure 2, define the ratio of transmitted diffuse intensity to the 
specular incident intensity and the ratio of the reflected diffuse intensity to the 
specular incident intensity. 


The haze functions are calculated by the expressions: 


4no-CT~-|n,—n|\"" 
H,; =1— exp -(——)} (1) 
4no ~ CT -n,\"* 
H, —=1— exp (= ) (2) 


Table 1. Default parameter values for the haze functions 


No. Parameter Value Units 
1 CR 1 
2 CT 0.5 
3 NR 2 
4 NT 3 
5 oO 20 nm 


where A is the optical wavelength, CT, CR, NT, NR and o are user definable 
model parameters. The parameter o is a measure of the mean feature size at the 
interface. The default values of the model parameters are given in Table 1. 


Figure 3 shows example calculations of the haze functions versus wavelength for 
various values of o. 


The diffuse light is reflected and transmitted in random directions from the 
interface. The distribution of diffuse light can be described by an angular 
distribution function, ADF. The ATLAS-Luminous simulator supports several 
user selectable ADFs as summarized in Table 2. 
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Fig. 3. Example haze functions as a function of wavelength for various 6 [7]. 


Table 2. Un-normalized Angular Distribution Functions in Luminous. 
No. ADF(®) 
1 Constant 1 
20 
9; Triangle i 
ir 
_¢* } 
3 Gauss exe spaeaar, 
aa 
4 Lorenz &XP\ e?, orspeasron® 
5 Lambert cos(®) 
cos(P) 
6 Ellipse a 0.5 pe 
cos*(®) + (semmmvon) sin?) 


In Table 2, ® is the angle of propagation with respect to the normal to the 
interface and the ADFs are un-normalized. Normalization insures conservation of 
intensity derived from the haze functions above. 


The geometrical contribution of the diffuse light is accounted for in the calculated 
photogeneration rate given in equation 3 where P, is the intensity at the interface, 
y is the distance normal to the interface and o is the absorption coefficient. 


2 sy 
Gly) =2 | P, ADF(®)aecostd® 
0 


(3) 


Reflection and refraction of diffuse light at subsequent interfaces is modeled using 
the Fresnel equations. 
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3. The simulation software — Silvaco ATLAS 


Atlas is a 2D and 3D device simulator that performs DC, AC, and transient 
analysis for silicon, binary, ternary, and quaternary material-based devices. Atlas 
enables the characterization and optimization of semiconductor devices for a wide 
range of technologies. 


Luminous 2D/3D is an advanced device simulator specially designed to model 
light absorption and photogeneration in planar and non-planar semiconductor 
devices. Solutions for general optical sources are obtained using geometric ray 
tracing or beam propagation methods. These features enable Luminous 2D/3D to 
account for arbitrary topologies, internal and external reflections and refractions, 
polarization dependencies, dispersion and coherence effects. Luminous 2D/3D is 
fully integrated within Atlas with a seamless link to other Atlas device technology 
modules [19,24]. 


Solar cell characteristics, such as collection efficiency, spectral response, open 
circuit voltage and short circuit current can be extracted with Luminous 2D/3D. 


4. Results and discussions 


A thin film hydrogenated amorphous silicon cell was simulated using the Silvaco 
ATLAS simulation software. The device consists of a one half micron thick 
textured transparent conducting oxide, on top of p/i/n textured layers of 
amorphous silicon of thicknesses 9nm/650nm/20nm all on top of 300 nm of 
silver, shown in Figure 4. 
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Fig. 4. The structure of a thin film a-Si:H solar cell 
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In the following, the steps of the simulation in ATLAS, together with the obtained 
results, are presented. 


First we specify a mesh of 0.1 square meter. Since the structure is essentially one 
dimensional we ignore variations if x and z. 

Next, we introduce a set of material layers comprising the solar cell device. It is 
important to introduce all material layers even if they have little or no electrical 
effect since they will have significant effects on the transfer matrix analysis used 
to model the optical problem. 

Next, we define several of the interfaces as diffusive with certain haze function 
parameter specifications and angular distribution functions. 


Also, we define the complex index of refraction for some of the active layers 
using tabular data. 

For the metals and oxide we use the Sopra data base for the default complex index 
of refraction [20, 21]. 


We plot the complex index of refraction data (see Figure 5). 
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Fig. 5. The plotted complex index of refraction data: 

a) real; b) imaginary. 
We model the optical source as coming from above the device normal to the top 
surface. We define the source as monochromatic with an arbitrary wavelength 
since we will be performing spectral analysis later. We also enable the transfer 
matrix with diffuse interfaces. 


We save the spectral response for the device to a log file. We can extract useful 
performance measures using the EXTRACT statements. 


Finally, we can plot the results as compared against measurement. Figure 6 shows 
excellent comparison between experiment and simulation for this device [22-24]. 
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(cathode current)/(source photo current) vs optical wavelength 
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Fig. 6. Quantum Efficiency comparison between experiment and 
simulation of a thin film a-Si:H solar cell. 


5. Conclusions 


This simulation demonstrates the use of the transfer matrix method, with diffusive 
interfaces to model a textured a-Si solar cell. This approach is equally applicable to other 
types of textured thin film cells such as CIGS cells. The keys to this simulation are the 
use of the transfer matrix method to account for the coherent interference in the thin films 
and the use of diffusive interfaces to account for the effects of light trapping in the 
textured surfaces. 


This example also shows the importance and methodology of including conducting layers 
in the simulation using the transfer matrix method. We have implemented a self- 
consistent optoelectronic model for simulation of solar cells with rough textured 
interfaces. 


This simulator combines simulation of specular and diffuse light with haze functions 
describing the relation between the light intensities and angular distribution functions 
describing the diffuse dispersion of light. 


The comparison with experiment gives good accuracy. 
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